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An exact Craik-Criminale (CC) solution to the incompressible Navier-Stokes (NS) equa- 
tions describes the instability of an elliptical columnar flow interacting with a single 
Kelvin wave. These CC solutions are extended to allow multi-harmonic Kelvin waves to 
interact with any exact "base" solution of the NS equations. The interaction is evaluated 
along an arbitrarily chosen flowline of the base solution, so exact nonlinear instability in 
this context is locally convective, rather than absolute. Furthermore, an iterative method 
called "WKB-bootstrapping" is introduced which successively adds Kelvin waves with 
incommensurate phases to the extended CC solutions. This is illustrated by constructing 
an extended CC solution consisting of several Kelvin waves with incommensurate phases 
interacting with an elliptical columnar flow. 



1. Introduction 

ICraik fc Criminalel l)l986|) showed formally that the sum of a base flow that is linear 
in the spatial coordinates and a single traveling wave is a global solution of the in- 
compressible Navier-Stokes (NS) equations on an unbounded spa tial domain. Although 
this fact was known in certain specific cases (see, for example, ( Chandrasekharlll96ll 



pp. 85-86)), it was Craik & Criminale who formalized the general framework. This dis- 
covery provided an exact nonlinear interpretation of elliptic instability, the resonant 
mechanism by which v ortex stretch ing creates three-dimensional instabilities in swirling 
two-dimensional flows <|Bavlvlll98^ . This nonlinear interpretation of elliptic inst ability 
also holds when additional physics such as buoyanc y IjMivazaki fc Fukumotdll992IL r ota- 
tion <|Craiklll989HMivazakilll99ilBavlv el oi.lll99^ and magnetic fields <|CrailJll988() are 
included. Mos t recently, elliptic in s tability has been used to st udy the effects of turbu- 
lence models (jCambon et aZ.lfl994HFabiionas fc Holmll2003Ql6T l. Elliptic instability and 
its interpretati on using the Craik-Criminale (CC) solutions were recently reviewed by 



its mtcrpretatic 
iKerswelll (120021) 



A second equally ground-breaking body of work related to the CC class of solutions 
is Lifschitz's WKB-like stability meth od for fluid mechanics known as the Geq metri- 
cal Optics Stability Method (GOSM) flLifschitz fc Hameiril Il99lt lLifschitzlll994[K The 
method linearizes the perturbation equations about the unperturbed flow and exam- 
ines the evolution of a high-frequency infinitesimal wavepacket in a Lagrangian reference 
frame which moves with the unperturbed flow. The power of this method is its insensi- 
tivity to the unperturbed flow's global geometry and time dependence. Remarkably, the 
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equations governing the nonzero leading order terms of GOSM coincide with those for 
CC solutions. This coincidence motivates the present paper. 

This paper bridges the gap between GOSM and CC solutions by extending the CC 
class of exact nonlinear solutions of the NS equations in two stages. First, we show that 
the CC construction method works when a Kelvin wave together with any number of 
its harmonics interacts with any exact solution of the nonlinear NS equations, provided 
the interaction is evaluated along some arbitrarily chosen flowline of the base solution. 
This extension of the CC construction met hod changes its interpretati on from absolute 
instability to locally convective instability ijHuerre fc Monkewit dll990|) . Second, we in- 
troduce an iterative method called "WKB-bootstrapping" which constructs extended CC 
solutions interacting with additional traveling waves whose phases are incommensurate. 
These are multi-frequency CC solutions. (The phrase "incommensurate phases" will be 
explicitly defined in Eq. I|3.2|l .) The WKB-bootstrapping algorithm allows Kelvin waves 
with incommensurate phases to be added to the base flow one at a time. The power 
of this method is illustrated by deriving the exact solution for nonlinear instability of 
an elliptical columnar flow that is interacting with several traveling Kelvin waves whose 
phases are incommensurate. Also, a steady circular columnar flow interacting with a sin- 
gle Kelvin wave is found to be critically stable. That is, either perturbing the streamline 
to introduce eccentricity, or adding multiple traveling waves with incommensurate phases 
will yield Kelvin waves with exponentially growing amplitudes. 



2. Craik-Criminale solutions for generalized base flows 

We consider the NS equations in a coordinate system rotating about the z— axis with 
angular velocity 17 = Cle z : 

d t u + u • Vu + 2Q. x u + Vp = fAu. (2.1) 

Here Vu = 0, where u is a velocity field, p is the pressure, and v is the kinematic viscosity. 
Let U be any exact solution to Eq. I|2.1|l with corresponding pressure P. Consider a new 
exact solution of the form u = U + U', p = P + P'. In the analysis of nonlinear instability 
to follow, we shall refer to U as the 'base' flow and U' as the 'disturbance.' We proceed 
by inserting this additive velocity decomposition into Eq. (|2.1|l . and using the condition 
that U is an exact solution, together with the vector identity 

(u- V)v = (V x v) x u+ (Vv) T • u. (2.2) 

Consequently, the nonlinear disturbance dynamics is governed by 

d t V + (U • V)U' + (VU) T • U' + (U' • V)U' + n x U' + VP' = vAXJ', (2.3) 

and V • U' = 0. Here, II = 2fi + V xU denotes the total vorticity of the base flow. We 
choose the disturbance to consist of a traveling wave and its harmonics of the form 

oo oc 

U'= M H a m (*)e im/3 *, P'= V |m| Pm(t)e im ^. (2.4) 

m— — oc m— — oo 

Here, $ = $(i, x) is the fundamental phase of the disturbance. In linear analysis, the sum 
U + U' would be a solution to the NS equations linearized about U; that is, Eq. (|2.3|l with 
the term quadratic in U' removed. Craik & Criminalc showed that if VU is a function 
of time only, then the sum U + U' is an exact solution to the nonlinear equations, where 
$ is chosen to be linear in x. We now extend the CC solutions, by showing the sum 
is an exact solution to the NS equations for any base flow U, when evaluated on an 
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arbitrarily chosen flowline of U. Specifically, we consider the flowline x = \{t, xo) where 
dx/dt = U(t, %) with x(0, xo) = xo, and we label this trajectory Q Xo . We emphasize that 
the choice of the flowline is arbitrary. However, the subsequent equations and analysis 
will be evaluated on this flowline. We seek solutions of Eq. I|2.3|) in the above forms whose 
phase $ is frozen into the base flow. That is, 



(d t + U • V) $ = 0, 



(2.5) 



with x) = k(t) • x + 8(t). Along <2 Xo , the frozen phase relation in Eq. (|2.5() has the 
solution $(t, x(t, xo)) = $(0, xo) = $o(xo)- The parameters fi, (3 are scaling factors which 
allow us to choose |k(0)| = 1 and max m {|a m (0)|} = 1. Evaluating the incomprcssibility 
condition V • U' = on Q Xo immediately yields 



V-U' 



J2 fi lml imp\a m (t)-k(t) 



0. 



Since the bracketed expression depends only on time along Q Xo , we conclude that 

a m • k = (2.6) 
for all m. Evaluating the quadratic term U' • VU' in Eq. I|2.3[) along Q Xo yields 



U' • VU' 



E 



m| 



a m (t) • k(i) 



J2 »n^ |n| a n (t)e in ' w 



which also vanishes exactly. Thus, on Q Xo , the incomprcssibility condition V ■ U' = 
implies the transversality condition a TO • k = 0. This, in turn, implies that the quadratic 
term U' • VU' in Eq. (|2.3() exactly vanishes on Q Xf) . Substituting the frozen-phase solution 
ansatz in Eq. (|2.5|) into Eq. (|2.3|) yields the Eulerian equation, 



(d t + (VU) T ) a m + n x a m - m/3p m V$ 



-z/TO 2 /3 2 |V$| 2 a„ 



(2.7) 



At this point, it is clear that we can take a m and p m to be real valued functions without 
loss of generality. Furthermore, it follows that both the real and imaginary parts of 
Eq. I|2.4|) individually are solutions. We may further decompose the base flow as U(x, t) = 
U(x,<) + C(t), where U(x ,i) = C(t), and set d t S + C • k = 0, where (5(0) = $ (x ) - 
k(0) • Xo. This relation eliminates all contributions of C(t) to the evolution of <& and 
uniquely implies 8(t). Thus, along Q Xo , Eq. I|2.3|l reduces to a set of ODEs. These are 
obtained by taking the gradient of Eq. (|2.5|) , and evaluating Eq. I|2.7|l along <2 X0 : 



dtX 

(dt- 

(d t - 

Here, the quantity 



= U, 

C T ) k = , 



C 1 a„ 



II x a rn — m/3p m k = — vm 2 f3 2 |k| 2 a„ 



Ci 



(VU), 



x=x(t,x ) 



diUi 



(2.8) 
(2.9) 
(2.10) 

(2.11) 



x=x(t,x ) 

is the velocity gradient tensor along Xi an d dt — dt + U ■ V is the material derivative. 
Note that, C = VU is a function of time only, when evaluated along Q Xo . The operator 
dt + C T comprises the total time deri vative in a Lagrang i an fra me that moves with the 



base flow U along Q Xo . In contrast. ICraik fc Criminald (jf 986|) considered the case of 
global Eulerian solutions, rather than solutions along particular Lagrangian trajectories. 
In their framework, only base flows of the form U = C(t)x + C(t) are admissible. We 
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have shown, however, that along a single, arbitrarily chosen Lagrangian trajectory, the 
CC construction method works for any exact solution of the NS equations. Our result 
reduces to Craik & Criminale's in the case when global Eulerian solutions are required. 

We rewrite the pressure p m in Eq. Ij2.1t)|l as follows. Upon taking the dot product of 
Eq. I|2.10p with k and noting that a m • k is an integral of motion, we solve for the pressure 
contribution p m of the m-th harmonic in terms of the amplitude a m and the wavevector 
k, as 

mf3p m \k.\ 2 = {£. + C T )a m ■ k + II x a ln • k. (2.12) 

Finally, the viscosity may be removed from the right hand side of Eq. (|2.10() by the 
change of variables, 

a m = a m cxp (- J vm 2 (3 2 \k(r)\ 2 dr^j . (2.13) 

One then solves for Eq. 12.1L)fl with the right hand side replaced by the zero 

vector. 

We emphasize that the construction method presented here is based on three assump- 
tions: (1) the disturbance amplitudes a m and p m are independent of the spatial variable; 
(2) 4> is linear in x = %(i,Xo) along a flowline in the Lagrangian frame of the base flow 
U; and (3) all of the analysis is carried out along a single, arbitrarily chosen flowline 
specified by its initial position xo. The construction method presented here fails if a m 
is allowed to have spatial dependence. Such spatial dependence may be admitted by 
linearizing the equations of motion about U and introducing high frequency wavepack- 
ets rather than trad itional traveling waves ( Lifschitz fc Hameirilll99ll iLifschitd Il99l 



lOiaia fc Hohnlll99fil) . We note that the superposition of traveling wave harmonics does 
not alter the evolution of the base flow. This is because the base flow U is an exact 
solution, by itself, whose evolution drives the disturbance in Eq. (|2.3(l . 



3. WKB bootstrapping 

The general framework outlined in the previous section allows for the construction of 
multi- frequency CC solutions. We call this method WKB-bootstrapping' and outline an 
iterative algorithm for applying it. Let us denote iterations of the WKB-bootstrapping 
method with a parenthetical superscript. Beginning with an exact base flow U = u^ ' 
for n = 0, we iteratively construct 

u (n+i) = u {n hJ2^ {n+1) y n ^ +1) sin(m/3 ( ™ +1) $ (n+1) ), (3.1) 

m 

where $( rl+1 ) = $( n+1 )(i,x), by following the process described above in Section 2. 
Specifically, at iteration n + 1, we take U = u(") and U' is the new traveling wave. 
The evolution of the solution at iteration n + 1 is carried out along a single flowline 
in a Lagrangian frame which moves with . That is, the phase at iteration n + 1 is 
(t, x) = k(" +1 ) (t) ■ x + (t), where x = x {n+1) (*, *o) with x (n+1) (0, x ) = x 

and dx {n+1) /dt = u( n \t,x {n+1) )- Thus, the velocity gradient tensor C {n+1) = Vu'">, 
evaluated at x = x' ri+1 - ) (t, xq), is again a function of time only, along the trajectory 
Qx™ ■ Though at each iteration, the trajectory Qx™ w i u differ from trajectories of 
previous iterations, it is crucial that all iterations pass through the same point Xo, that 
is> (0, xo) = xo, j = 1, 2, . . . , n + 1. The construction method allows the phase $( n+1 ) 
of the new Kelvin wave to be incommensurate with the previous phases, by which we 
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mean that the phases are not functionally related for all x and t. That is to say, 

$(" +1 )(t,x)^/(<f>«(t,x)), j = l,2,...,n (3.2) 

for all x, t and any function /. Alternatively, the above definition for incommensurate 
phases can be rewritten as 

k (" +1 )(i)^/'($«)k«(t), i = l,2,...,n (3.3) 

for all time. That is to say, none of the wave vectors are always parallel, although some 
wave vectors may be parallel at certain instants in time. 

Because and all of the evolve along different trajectories, the phases will 

tend to be incommensurate. This iterative process may continue ad infinitum. We em- 
phasize that the traveling wave at iteration n + 1 is affected by all the preceding waves, 
that is, those from iterations 1,2, ... ,n. However, the wave at iteration n + 1 does not 
influence the previous waves. 

We emphasize that it is crucial that the waves in the WKB-bootstrapping algorithm 
be added one at a time. For example, the solution u*- 3 -* = u' ' + /i^a^ sin^ 1 )^ 1 )) + 
^(2) a (2) s in(/3( 2 )$( 2 )) + ^( 3 )a( 3 ) sin(/3( 3 )<I>( 3 )) is not in general an exact global solution of 
the NS equations. In particular, one cannot conclude that Eq. (|2.(i|) holds for the waves 
individually. The construction method outlined in Section 2 relies on the fact that we 
build the solution in three steps. First, we construct = u(°) + /x^a^ sin^ 1 )^ 1 )) 
by setting U = and U' = fj^a' 1 ) sin(/3W$W) into the theory of Section 2. We 
are now confident that u^ 1 ^ is an exact solution in a Lagrangian frame of an arbitrar- 
ily chosen flowlinc of u^ ) which passes through some point Xo. Next, we construct 
u (2) = u (0) + M (i) a ( 1 )sin(/?( 1 )$( 1 )) + ^a^sin^ 2 )^ 2 )) by setting U = u« and 
U' = n^a!^ sin(/3( 2 )$ ( ' 2 ' ) ) into the same theory. Since U is an exact solution, it is 
removed from the problem for the evolution of a^ 2 ) and k^ 2 ). Then, the desired solu- 
tion is an exact solution along the flowlinc in the Lagrangian frame of u' 1 ) which passes 
through Xo. Finally, we construct the desired solution u^ 3 ) by setting U = u' 2 ) and 
U' = pt( 3 )a< 3 ) sin(/3( 3 )$( 3 )) into the same theory. This final solution is an exact solution 
in a Lagrangian frame of the flowline of which passes through xq. The individual 
flowlines are illustrated in Fig.^ Finally, we note that by adding the waves one at a time 
in different Lagrangian frames, we generate the transversality conditions a' 1 ' • k' 1 ) = 0, 
a ( 2 ) .k^ 2 ) = 0, and a^ 3 ) -k^ = needed in the CC construction method. These sequential 
incomprcssibility relations enable the successful application of WKB-bootstrapping, by 
eliminating the quadratically nonlinear term in the evolution equation for each frequency. 



4. Example: stability of a circular columnar flow 

We illustrate the WKB-bootstrapping method with an example. 

4.1. Primary instability 

Consider a base flow of the form u^ ^ = 5x whose streamlines are ellipses with eccentricity 
7 in a non-rotating coordinate system (fl = 0): 

/ -1 + 7 0\ 
5= 1 + 7 . (4.1) 
\ 0/ 

We first construct a classical CC flow of the form u' 1 ' = 5x + ^ 1 'a' 1 ' sin(/3W This 
corres ponds to the classic problem of elli ptic instability which was investigated bv lKelvinl 
j 18871) for the circular case (7 = 0) and bv lBavlvl <Jl98fil) for the elliptic case (7 > 0) almost 
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Figure 1. Representative graphic of the WKB-bootstrapping algorithm. First, 
xS 1 ' = u l °' + fv- 1 ' aS 1 ' sin(/3 ( ' 1 ^$* 1 ' ) ) is an exact solution in the Lagrangian frame 
with moves along an arbitrarily chosen flowline of u' ' (left figure). Second, 
u (2) _ u (0) + ^(i) a (i) sin(/3 (1) $ (1) ) + fi {2) SL {2) sin(/3 (2) <I> (2) ) is an exact solution in a Lagrangian 
frame which moves along a flowline of u' 1 ' which passes through xq (solid line in center figure). 
Finally, = + sin(/3«i>«) + M < 2 >a< 2 ) sm{p^^) + u^r^ sin(f3^^) is an 

exact solution in a Lagrangian frame which moves along a flowline of u' 2 - 1 which passes through 
xo (solid line in the right figure). 



a century later. The equation for the wave vector has the analytical solution = 
[sin 9 cost, ks>vh9 smt, cos 6] T , where 8 is polar angle the wave vector makes with the axis 
of rotation, n 2 = (1 — 7)/(l + 7), and i = t\/l — 7 2 . The equation for the amplitud e 
a' 1 -* with = S in turn satisfies a Floquet problem l|Yakubovich fc Starzhinskiilll967lb 
That is, E q. (|2~TT)jl can be written as d t a {1) = 7V(i)a (1) where Nit + r) = N(t) and 
t = 2,7r / y/l — 7 2 . Thus, Lyapunov-like growth rates are determined by computing the 
monodromy matrix V , which is the fundamental solution matrix with identity initial 
condition evaluated at t = r. The growth rat es are given by a = lnjmaxi |5R(pi)|}/r, 
where pi,i — 1,2,3 arc the eigenvalues of V . iBavlvl l|l986|) was able to show that for 
certain values of 7 and 6, the amplitude could grow exponentially in time. See the left 
figure in Fig . The linearized perturbation analysis (about u^ ) = Sx) performed by 
Bavlvl lll986l) concluded t hat the base flow u(°) was unstable for all 7 > 0. The work of 



Craik fc CriminaTel l)l986|) showed that Bayly's analysis was an exact solution. 



4.2. Secondary instability 

We now construct the next WKB-bootstrap iteration as vv ' = Sx+ziWaW sin(/3( 1 )$( 1 )) + 
^( 2 ) a ( 2 ) sin(/3( 2 )$( 2 )). The equations for the phase $( 2 ' and the amplitude a^ 2 ) are given in 
Eqs. 1)2.9(1 - 1(2.12(1 with U = u' 1 '. In the interest of simplicity, we consider the case cos£? = 
0.2,7 = 0. Then u' 1 ) is periodic with period 2tt, and the equation for a^ 2 ) again satisfies 
a Floquet problem. We further simplify our analysis by examining the flow at the origin, 
which is a fixed point of the system. Then, $W = and £ (2 *> = S + /wW/J^a^ (kW) T . 
Note that although $W = and <&( 2 ) = at this fixed point, this does not imply that the 
phases are commensurate. Indeed, we choose the initial condition for the wave vector 

k (2) 

to be arbitrary in three-dimensions: k^ 2 ^(0) = [cos ^ sin 0, sin^sin0,cos(/)] T . The center 
figure in Fig. [5] shows the domain in the (cos <fi, cos ip) plane for which the amplitu de 
of a' 2 ' grows exponentially. This figure first appeared in lLifschitz fc Fabiionasl l|l996l) in 
which the authors examined the st ability of Bayly's solution for 7 = using GOSM. We 
conclude, then that the analysis of lLifschitz fc Fabiionasl l)l996j) is also an exact solution 
of the full nonlinear Euler equations. This result applies to the subsequent articles by 
the same authors and their collaborators <)Fabiionas et all Il997t iFabiionas fc Lifschitzl 
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Figure 2. Left most figure: Initial condition parameter plane for u' 1 . For 7 = 0, the Kelvin 
wave is periodic in time; for values of 7, cos 9 which fall into the white region for 7 > 0, the 
Kelvin wave is quasiperiodic, that is, the periods of a' 1 ) and are rationally related; for 
the remaining region, the Kelvin wave has an exponentially growing amplitude. Center figure: 
Initial orientation parameter plane for u' 2 ' at x = with parameter values cos 9 = 0.2, 7 = 
and = 1. For parameter values which fall into the white regions, the amplitude of 

a' 2 ' is bounded; for all others, the amplitude of a' 2 ' grows exponentially in time. Right figure: 
Initial orientation parameter plane for u' 3 -* at x = with parameter values cos 9 = 0.2, 7 = 0, 
cosci = 0.5, cos ip = 0. For parameter values which fall into the white regions, the amplitude of 
a' 3 - 1 is bounded; for all others, the amplitude of aS 3 ^ grows exponentially in time. 



Il998(h The authors viewed the perturbation analysis of u' 1 ' as a secondary perturbation 
analysis of u' ' . 

4.3. Tertiary instability 

We construct a third iteration of this method as u^ 3 ) = 5x + /iWa' 1 ' sin(/9W$W) + 
^(2) a (2) s in( / g( 2 )<i)( 2 ))-|-^( 3 )a( 3 ) sin(/?( 3 )$( 3 )), where the phases are incommensurate. Again, 
we examine the case cos# = 0.2, 7 = at the fixed point x = 0, and additionally 
we choose cose!) = 0.5 and cost/i = 0. Thus, the flow u^ 2 ) is bounded in time, and 
£^ = S + /iWjjWaW^ 1 )) 7 + /^ 2 )/?( 2 )a( 2 )(k( 2 )) T The individual components of u< 2 ) 
are periodic, but there is no apparent common period. To analyze the evol ution of a^ 3 ) , 
we use a quasi-periodic analogue of Floquet theory as discussed in detail in iBavlv et all 
l)l996() . We use the incompressibility condition a^ 3 ^ • k^ 3 ) = to reduce Eq. (|2.10|l for a^ 3 ) 
to a two-dimensional system. We then use the Priifer transformation 

(i) = e d « cos[c(t)] , (i) = e d{t) sin[d(t)] , (4.2) 

where and a^ 3 ' are the first two components of a^ 3 ) . This produces a coupled pair of 
ordinary differential equations for c(t) and d(t). Note that linear growth of d(t) implies 
exponential growth for a' 3 ' . One can define two quantities called the growth rate / and 
the winding number W as 

/ = lim [d(t)/t], W= ]hn[c(t)/t], (4.3) 

t — >oc t — >oo 

respectively. It was first shown conclusively in IBavlv et al\ l)l996(l that these quan- 
tities my be approximated by long time simulations independently of the values of 
c(0) and d(0). We choose the initial condition for k 1 - 3 ) in polar coordinates as before: 
k( 3 )(0) = [cos x sin 77, sin % sin 77, cos 77] T . The right figure in Fig. |3 shows the parameter 
plane (cos rj, cos x) f° r which the a mplitude a^ 3 ) grows exponenti ally in time for a specific 
orientation of a^. In the spirit of iLifschitz fc Fabiionasl l)l996(l . we may view this as a 
tertiary perturbation of the rigidly rotating flow which again provides an exact nonlinear 
NS solution. 
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5. Discussion 

The present analysis extends the pre vious results f or stability of a rigidly rotating col- 



stable in the presence of a traveling wave, slight perturbations in the streamline eccen- 
tricity yield a traveli ng wave whose amplitude gro ws exponentially in time (left figure in 
Fig.0). Furthermore. lLifschitz fc Fabiionasl l)l996|) showed that while the primary travel- 
ing wave was always stable, the amplitude of a secondary wave with an incommensurate 
phase grows exponentially (center figure in Fig- El - In the present paper, we show that 
even when the secondary wave is stable (that is, then the amplitude is bounded for all 
time), the tertiary wave is unstable for a larger portion of its parameter space. Thus, 
we conclude that a circular columnar flow is critically stable in the sense that either 
perturbations in the eccentricity or the addition of incommensurate phases will result in 
waves whose amplitudes grow exponentially in time. 

This paper has accomplished t hree main objecti ves. First, we extended the pioneering 
work o flCraik fc Criminalel ilQSfift and lBavlvl l|l98fil) to allow nonlinear instability analysis 
of any exact solution of the NS equations due to multi-harmonic Kelvin waves, evaluated 
along a flowlinc in the Lagrangian frame of the base flow. Second, we showed how to 
construct a sequence of exact solutions to the NS equations that involves the multi- 
frequency superposition of Kelvin traveling waves by adding the waves one at a time. This 
iterative method of constructing exact NS solutions was called 'WKB-bootstrapping.' 
Third, we used the WKB-bootstrapping method to unite and extend various classical 
linearized instability analyses of a circular columnar flow. 

The authors arc indebted to A. Lifschitz-Lipton whose work has guided us throughout 
this investigation. We are also grateful to anonymous referees whose comments signif- 
icantly improved the presentation of the material. BF thanks the Theoretical Division 
at the Los Alamos National Laboratory for their hospitality and acknowledges valuable 
exchanges with members of the Laboratory's Turbulence Working Group. 
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